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ABSTRACT 

The  lecture  introduces  the  basic  principles  of  hydrodynamic  lubrication  and  the  fundamental  equation  of 
Classical  Lubrication  Theory.  The  analysis  proceeds  to  derive  the  static  and  dynamic  performance 
characteristics  of  short  length  cylindrical  journal  bearings,  with  application  to  the  dynamic  forced 
performance  of  a  rigid  rotor  supported  on  plain  bearings.  In  a  radial  bearing,  the  Sommerfeld  number 
defines  a  relationship  between  the  static  load  and  the  journal  eccentricity  within  the  bearing.  This  design 
parameter  shows  the  static  performance  of  the  bearing  as  rotor  speed  increases.  Rotordynamic  force 
coefficients  are  in  troduced  and  their  effect  on  the  stability  of  a  rotor-bearing  system  thoroughly  discussed. 
Cross-coupled  force  coefficients  are  solely  due  to  journal  rotation,  and  the  magnitude  ( and  sign )  of  the 
cross-stiffness  determines  rotordynamic  stability.  The  whirl  frequency  ratio  (WFR)  relates  the  whirl 
frequency  of  sub  synchronous  motion  to  a  threshold  speed  of  instability.  The  desired  WFR  is  null;  however, 
plain  cylindrical  bearings  show  a  whirl  ratio  of  just  0.50,  limiting  the  operation  of  rotating  machinery  to 
shaft  speeds  below  twice  the  system  first  critical  speed.  The  analysis  concludes  with  a  review  of  practical 
(in  use )  journal  bearing  configurations  with  highlights  on  their  major  advantages  and  disadvantages, 
including  remedies  to  reduce  or  entirely  avoid  subsynchronous  whirl  instability  problems. 

1.0  FUNDAMENTS  OF  FLUID  FILM  BEARING  ANALYSIS 

Figure  1  depicts  an  idealized  geometry  of  a  fluid  film  bearing.  The  major  characteristic  of  a  lubricant  film, 
and  which  allows  a  major  simplification  of  its  analysis,  is  that  the  thickness  of  the  film  (h)  is  very  small 
when  compared  to  its  length  (L)  or  to  its  radius  of  curvature  (R),  i.e.  (h/L)  or  (h/R)  <«  1. 


Figure  1 :  Geometry  of  Flow  Region  in  a  Thin  Fluid  Film  Bearing  (h  «  Lx,  Lz). 

San  Andres,  L.  (2006)  Hydrodynamic  Fluid  Film  Bearings  and  Their  Effect  on  the  Stability  of  Rotating  Machinery.  In  Design  and  Analysis 
of  High  Speed  Pumps  (pp.  10-1  -  10-36).  Educational  Notes  RTO-EN-AVT-143,  Paper  10.  Neuilly-sur-Seine,  France:  RTO.  Available 
from:  http://www.rto.nato.int/abstracts.asp. 


RTO-EN-AVT-143 


10-1 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1 .  REPORT  DATE  2.  REPORT  TYPE 

01  NOV  2006  N/A 

3.  DATES  COVERED 

4.  TITLE  AND  SUBTITLE 

Hydrodynamic  Fluid  Film  Bearings  and  Their  Effect  on  the  Stability  of 
Rotating  Machinery 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Mast-Childs  Tribology  Professor  Turbomachinery  Laboratory  Texas 

A&M  University  College  Station,  TX  77843-3123  USA 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

See  also  ADM002051.,  The  original  document  contains  color  images. 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

18.  NUMBER  19a.  NAME  OF 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  |J|J 

unclassified  unclassified  unclassified 

36 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Hydrodynamic  Fluid  Film  Bearings  and 

Their  Effect  on  the  Stability  of  Rotating  Machinery 


ORGANIZATION 


A  plain  cylindrical  journal  bearing,  see  Figure  2,  comprises  of  an  inner  rotating  cylinder  (journal)  of 
radius  Rj  and  an  outer  cylinder  (bearing)  of  radius  RB  (>Rj).  The  two  cylinders  are  closely  spaced  and  the 
annular  gap  between  the  two  cylinders  is  filled  with  some  lubricant.  The  radial  clearance  c  -  (RB  -Rj)  is 
very  small.  In  most  fluid  film  bearings  with  incompressible  liquids,  c/RB  =0.001;  while  for  gas  film 
bearings,  C/RB  = 0.0001 ,  typically. 


As  a  consequence  of  the  smallness  in  film  thickness,  the  effects  of  the  film  curvature  are  negligible  in 
the  operation  of  a  journal  bearing. 

The  analysis  of  the  flow  equation  defines  the  circumferential  flow  Reynolds  number  (Re)  as 


„  pU *  c 

Re  =  — -  &  Rc  =Rc  (c/L*) 

A 


(la) 


based  on  the  characteristic  speed  ( U*=QR ).  Re  denotes  the  ratio  between  fluid  inertia  (advection)  forces 
and  viscous-shear'  forces.  And, 


Re,  = 


pcoc~ 

A 


(lb) 


is  the  squeeze  film  Reynolds  number  representing  the  ratio  between  temporal  fluid  inertia  forces  due  to 
transient  motions  at  a  characteristics  frequency  (co)  and  viscous-shear  forces.  Fluid  inertia  effects  in  thin 
film  flows  are  of  importance  only  in  those  applications  where  both  Reynolds  numbers  are  larger  than 
ONE1,  i.e.  Re*,  Res  >>  1. 


The  film  thickness  to  characteristic  length  ratio  (c/L*)  in  thin  film  flows  is  typically  very  small.  Thus,  fluid 
inertia  terms  are  to  be  retained  for  flows  with  Reynolds  numbers  of  the  order  [2] : 

Re  >  (L*/c)  x  1  >  1,000  for  (c/L*)  =  0.001  (2) 

Classical  lubrication  is  based  on  the  assumption  that  fluid  inertia  effects  are  negligible,  i.e.  (Res,  Re*)  — >  0, 
in  most  practical  applications;  and  hence,  rendering  effectively  an  inertialess  fluid. 


In  actuality,  Re  >  12  for  steady  super  laminar  flow  in  thin  film  bearings,  as  demonstrated  in  [1].  In  thin  film  flows,  transition  to  turbulence  is 
due  to  instability  of  shear  driven  parallel  flow.  A  transition  to  turbulence  initiated  by  the  appearance  of  Taylor  vortices  generated  by 
centrifugal  forces  is  more  peculiar  to  configurations  with  large  clearances,  i.e.  not  common  in  thin  film  bearings.  The  accepted  Reynolds 
number  for  flow  turbulence  in  journal  bearings  is  R ec=pQRc//j  >2,000  [2]. 
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The  thin  film  laminar  flow  of  an  incompressible,  inertialess  and  isoviscous  fluid  is  governed  by  the 
following  equations  continuity  and  momentum  transport  equations: 


3(0  !  d(v-y)  {  d(vJ^Q 

8x  8y  dz 


5  P 

8  x 


+  jU 


0  = 


8P 

8  z 


+  M 


(3) 


where  vxyiZ  are  the  velocity  components  and  the  film  pressure  (P)  is  uniform  across  the  film  thickness  (/;). 
The  momentum  equations  establish  a  quasi-static  balance  of  pressure  and  viscous  forces.  The  lubrication 
equations  denote  “slow”  flow  conditions  with  time  appearing  as  a  parameter,  not  an  independent  variable. 

Table  1  presents  the  circumferential  flow  Reynolds  numbers  (Re)  for  a  typical  journal  bearing  application 
operating  with  different  fluids.  The  example  bearing  is  a  3  inch  (76  mm)  diameter  ( 2Rj )  journal  and  the 
clearance  to  radius  ratio  ( c/Rj )  is  0.001,  a  typical  value  for  journal  bearings.  Two  rotational  speeds  of 
1,000  and  10,000  lpm  (£2=104.7  and  1047  rad/s)  are  noted  in  the  table.  The  calculated  Reynolds  numbers 
(Re)  show  that  bearing  applications  with  mineral  oils  and  (even)  air  do  not  need  to  include  fluid  inertia 
effects,  i.e.  Re  <  1,000.  However,  process  fluid  applications  using  water,  R134a  refrigerant  and  cryogenic 
fluids  show  large  Reynolds  numbers  at  a  speed  of  10,000  rpm. 


Table  1 :  Importance  of  Fluid  Inertia  Effects  on  Several  Fluid  Film 
Bearing  Applications  (c/Rj  )=0.001,  Rj=  38.1  mm  (1.5  inch) 


Fluid 

Absolute 
viscosity  (ji) 
lbm.ft.s  x  10'5 

Kinematic 
viscosity  ( v) 
centistoke 

Re  at  1,000  rpm 

Re  at  10,000  rpm 

Air 

1.23 

15.4 

9.9 

99 

Thick  oil 

1,682 

30.0 

5.1 

51 

Light  oil 

120 

2.14 

71 

711 

Water 

64 

1.00 

159 

1,588 

Liquid  hydrogen 

1.075 

0.216 

705 

7,052 

Liquid  oxygen 

10.47 

0.191 

794 

7,942 

Liquid  nitrogen 

13.93 

0.179 

848 

8,477 

R134  refrigerant 

13.30 

0.163 

930 

9,296 

Note  that  current  bearing  applications  using  process  liquids  to  replace  mineral  oils  may  operate  at  speeds 
well  above  10,000  rpm.  Incidentally,  the  operating  speed  of  cryogenic  turbopumps  is  on  the  order  of  30  - 
70  krpm,  and  future  applications  (currently  in  the  works)  will  operate  at  speeds  close  to  200  krpm! 

Incidentally,  process  gas  and  liquid  seals,  isolating  regions  of  high  and  low  pressures  in  a  typical 
compressor  or  pump,  have  larger  radial  clearances  than  load  support  fluid  film  bearings.  For  example,  in 
water  neck-ring  and  interstage  seals  in  pumps,  RJC  ~  250,  and  thus  fluid  inertia  effects  are  of  importance 
even  at  relatively  low  rotational  speeds  (-1,000  rpm  and  larger).  The  topic  of  seals  is  analyzed  in  the  next 
lecture. 

1.1  Other  Fluid  Inertia  Effects 

Fluid  inertia  within  thin  film  flow  domains  can  be  safely  ignored  in  most  conventional  (oil  lubricated) 
bearing  applications.  However,  fluid  inertia  effects  may  also  be  of  great  importance  at  the  inlet  to  the  film 
and  discharge  from  the  film  sections  in  a  typical  pad  bearing  or  seal  application,  see  Figure  3.  Depending 
on  the  flow  conditions  upstream  of  a  sudden  contraction  or  a  sudden  enlargement,  a  fraction  of  the 
dynamic  pressure  head,  typically  given  as  (V2  p  U2 ),  is  lost  or  recovered. 
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U 

- ► 


Figure  3:  Pressure  Drop  &  Rise  at  Sudden  Changes  in  Film  Thickness. 


Sudden  pressure  losses  are  typical  at  the  edges  of  a  pocket  in  a  hydrostatic  bearing  and  at  the  inlet  plane  of 
annular  pressure  seals.  The  same  phenomenon  also  occurs  at  the  leading  edge  of  a  bearing  pad  in  high 
speed  tilting  pad  bearings.  A  sudden  pressure  recovery  is  also  quite  typical  at  the  discharge  section  of  a 
pressurized  annular  or  labyrinth  seal.  Note  that  the  importance  of  fluid  inertia  effects  may  be  restricted 
only  to  the  inlet  and  discharge  sections,  and  may  not  be  relevant  within  the  thin  film  flow  domain. 


2.0  REYNOLDS  EQUATION  AND  KINEMATICS  OF  JOURNAL  MOTION 

Lubricated  cylindrical  bearings  are  low  friction,  load  bearing  supports  in  rotating  machinery.  These  fluid 
film  bearings  also  introduce  viscous  damping  that  aids  in  reducing  the  amplitude  of  vibrations  in  operating 
machinery.  Figure  4  shows  a  schematic  view  of  a  cylindrical  bearing.  The  journal  spins  with  angular  speed 
(D.)  and  its  center  (Oj),  due  to  dynamic  loads,  also  describes  translational  motions  within  the  bearing 
clearance.  The  bearing  or  housing  is  stationary  in  most  applications.  Notable  exceptions  are  those  of 
floating  ring  journal  bearings  and  automotive  reciprocating  engine  support  rod  bearings. 
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x=R© 


Note:  film  gap  enlarged  for  description 
purposes 


Bearing 

0=9+([) 

Clearance  C=Rb-Rj 


Figure  4:  Schematic  View  of  a  Cylindrical  Journal  Bearing. 
Fixed  Coordinate  Systems  (0,z)  and  Moving  Coordinate  System  (0,z). 


2.1  Reynolds  Equation  for  Journal  Bearings 

The  smallness  of  this  ratio  allows  for  a  Cartesian  coordinate  (x=R0,  y,  z)  to  be  located  on  the  bearing 
surface  (see  Figure  4).  Then,  in  Classical  Lubrication,  Reynolds  equation  describes  the  generation  of 
hydrodynamic  pressure  (P)  within  the  bearing.  This  equation  arises  from  integration  of  the  momentum 
equations  (3)  across  the  film  thickness  and  substitution  into  the  continuity  equation  [2]: 


A{  ph}^{P 

8t  230  ^"30[  12 //3©j  3z[l2//3zj 


(4) 


in  the  flow  domain  {0  <  0  <2n ,  -Vi  L  <z  <Vi  Lj,  where  h( 0  z„  t)  is  the  film  thickness,  L  is  the  bearing 
axial  length,  U  =  fJRj  is  the  journal  surface  speed,  and  (p,  p)  denote  the  lubricant  density  and  viscosity, 
respectively. 
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The  boundary  conditions  for  the  hydrodynamic  pressure  in  a  plain  cylindrical  bearing  are2: 

a)  the  pressure  is  continuous  and  periodic  in  the  circumferential  direction,  i.e. 

P(0,Z,t)  =  P(&+27T,Z,t)  (5) 

b)  the  pressure  equals  the  discharge  or  atmospheric  value  (Pa)  on  the  bearing  sides,  i.e. 

P(  0,  Vi  L  ,  t)  =  P(  0,  -  Vi  L,  t)  =  Pa  (6) 

As  a  constraint,  the  hydrodynamic  pressure  needs  to  be  greater  than  the  liquid  cavitation  pressure 
everywhere  in  the  flow  domain,  i.e. 

P  >Pcav  in  0  <  0<2n,  -ViL  <z  <ViL  (7) 

Here  Pcav  represents  the  lubricant  saturation  pressure  or  the  ambient  pressure  needed  for  release  of 
dissolved  gases.  In  practice,  no  distinction  is  made  between  these  two  values  since  hydrodynamic  film 
pressures  could  be  one  to  two  orders  of  magnitude  larger  than  the  ambient  value. 

Consider  the  journal  and  bearing  to  be  aligned  and  the  journal  center  to  have  an  eccentricity  displacement 
e  (<c).  The  film  thickness  is 


h=c  +  ecos(0)  (8) 

This  formula  is  accurate  for  (c/7?)  ratios  as  large  as  0. 10.  The  film  thickness  derived  assumes  rigid  bearing 
and  journal  surfaces,  uniform  axial  and  azimuthal  clearance  and  no  journal  misalignment. 

2.2  Kinematics  of  Journal  Motion 

The  journal  center  Oj  is  displaced  a  distance  (e)  from  the  bearing  center  0B  ■  This  distance  is  known  as  the 
journal  eccentricity  and  may  vary  with  time  depending  upon  the  imposed  external  load  on  the  bearing. 
The  journal  eccentricity  cannot  exceed  the  bearing  clearance,  otherwise  solid  contact  and  potential 
catastrophic  failure  may  occur.  The  eccentricity  components  in  the  ( X ,  Y)  fixed  coordinate  system  are: 

ex-e  cos(  (j>);  eY-e  sin(  (j>)  (9) 

where  cf>  is  known  as  the  bearing  attitude  angle,  and  0=O+(j).  Then,  the  film  thickness  also  equals 

/?  =  c  +  ex  cos0  +  ey  sin©  =  e  sin#  (10) 

8  h  8  h 

and  - =  -ex  sin 0  +  eY  cos 0;  —  =  ex  cos0  +  ey  sin©  (11) 

00  8t 

where  ( .  )  denotes  differentiation  with  respect  to  time,  i.e.  (0/  8t). 

Substitution  of  the  film  thickness  gradients  into  Reynolds  equation  (4)  gives  the  following  lubrication 
equation  for  an  incompressible  and  isoviscous  fluid: 


The  simple  journal  bearing  model  does  not  account  for  feeding  holes  or  axial  grooves  for  supply  of  the  lubricant  into  the  bearing.  A  more 
detailed  discussion  on  lubricant  cavitation  and  its  physical  model  can  be  found  in  [3]. 
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1  8  j  h2  ap  j+  8  J  li2  dP\_{^ 


R2  8®  [12//  30 j  8z\\2ju8z 


nl  r. 

=  \ex  +eY  —  IcosO  +  <eY 


(12) 


An  alternative  form  of  Reynolds  equation  arises  when  using  the  angular  coordinate  (0).  This  angle  starts 
from  the  location  of  maximum  film  thickness.  A  coordinate  system  with  radial  and  tangential  (r,  t)  axes  is 
conveniently  defined  with  the  unit  radial  vector  along  the  line  joining  the  bearing  and  journal  centers. 

Recall  that  ex  -  e  cos( </>);  eY-e  sin( (/>),  and  e2  =  e\  +  eY . 
The  journal  center  velocities  in  the  (X,  Y)  and  (r,  t) 
coordinate  systems  are  related  by  the  transformation: 


cos  (j)  -sin  (f> 
sin  (j)  cos  (j) 


e 

e(j) 


(13) 


Note  that  Vr=e,  Vt=e</>  are  the  radial  and  tangential 
components  of  the  journal  translational  velocity.  From  the 
film  thickness  h  =  c  +  e  cos  0  ,  it  follows 


Translational  velocities  of  journal  center 


8h 

~8t 


—  =  -esin# 

86 


d6 

ecosO  -e  —  sin#-  =  ccos 9  +  ed)  sin#-  =  Vr  cos#  +  V,  sin# 
8t 


(14) 


Thus,  Reynolds  equation  (4)  for  an  incompressible  and  isoviscous  fluid  is  also  expressed  as 


1  8  f  h3  8P\  8  \  h 3  8P 


R2  86  {\2/u89\  8z\ll/u8 


=  ecos#  +e\<j)  _~r  sin^ 


(15) 


Equation  (15)  is  of  particular  importance  since  it  allows  us  to  realize  an  important  physical  phenomenon. 
Consider  the  journal  center  to  describe  circular  centered  orbits  with  a  fixed  amplitude  or  radius,  e.  Flence 
dc/dt=0.  Furthermore,  if  the  frequency  of  whirl  equals  to  50%  of  the  rotational  speed;  (j)  =  Q/2 ;  then  the 
right  hand  side  of  Eqn.  (15)  is  null;  and  hence  the  pressure  is  zero,  P=0  within  the  bearing  film  land. 
There  is  no  generation  of  hydrodynamic  pressure,  thus  resulting  in  a  sudden  loss  of  load  support 
capability.  The  50%  sped  whirl  phenomenon  is  the  basis  of  rotordynamic  instability,  as  explained  later. 


2.3  Bearing  Reaction  Forces 

Once  the  pressure  field  is  obtained,  fluid  film  forces  acting  on  the  journal  surface,  see  Figure  5,  are 
calculated  by  integration  of  the  pressure  field  acting  on  the  journal  surface.  An  equal  opposing  force  acts 
on  the  bearing  as  well.  The  bearing  reaction  forces  are  expressed  in  the  fixed  (X,  Y)  coordinate  system  and 
moving  (r,  t)  coordinate  system  as 


RTO-EN-AVT-143 


10-7 


Hydrodynamic  Fluid  Film  Bearings  and 

Their  Effect  on  the  Stability  of  Rotating  Machinery 


ORGANIZATION 


L  2tt  L  2  n 


Figure  5:  Fluid  Film  Force  Acting  on  Journal  Surface. 


The  relationship  between  the  fluid  film  forces  in  both  coordinate  systems  is  given  by: 


Fx 

Fy 


COS0 

-sin^ 

sin^ 

COS0 

F,_ 

(16. b) 


The  fluid  film  forces  are  generic  functions  of  the  journal  rotational  speed  (Cl)  and  the  journal  center 
translational  velocities,  i.e. 


Fa  =  Fa(Cl,ex,eY)=FQ 


f 

[  :  Q~P) 

e,e 

V 

2  Jy 

;  a  =  X,  Y  or  r,  t 


(17) 


An  analytical  solution  of  Reynolds  equation  for  arbitrary  geometry  cylindrical  bearings  is  not  feasible. 
Most  frequently,  numerical  methods  are  employed  to  solve  Reynolds  equation  and  to  obtain  the 
performance  characteristics  of  bearing  configurations  of  particular  interest. 

There  are  analytical  solutions  to  Reynolds  equation  applicable  to  two  limiting  geometries  of  journal 
bearings.  These  are  known  as  the  infinitely  long  and  infinitely  short  length  journal  bearing  models  [2], 

In  the  LONG  BEARING  MODEL,  see  Figure  6,  the  length  of  the  bearing  is  very  large,  L/D  —>  oo,  and 
consequently  the  axial  flow  is  effectively  very  small,  i.e.  ( dP/di )  =  0. 
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L 


Figure  6:  The  Long  Bearing  Model. 


For  large  L/D  ratios,  Reynolds  equation  reduces  to: 

i.W+£LL{*} 

R 2  d0{\2^  80}  8tl  2  30  1  J 


(18) 


This  bearing  model  gives  accurate  results  for  journal  bearings  with  slenderness  ratios  (L/D)  >  2.  Most 
modern  bearings  in  high  performance  turbomachinery  applications  have  a  small  L/D  ratio,  rarely 
exceeding  one.  Thus,  the  infinitely  long  journal  bearing  model  is  of  limited  current  interest.  Refer  to  [2] 
for  details  on  the  analytical  solution  of  Eqn.  (18). 

This  is  not  the  case  for  squeeze  film  dampers  (SFDs),  however,  since  the  long  bearing  model  provides  a 
very  good  approximation  for  tightly  sealed  dampers  even  for  small  L/D  ratios  [4] . 


3.0  STATIC  LOAD  PERFORMANCE  OF  SHORT  LENGTH  BEARINGS 

In  this  most  useful  bearing  model,  see  Figure  7,  the  bearing  length  is  short,  L/D— >  0,  and  consequently  the 
circumferential  flow  is  effectively  small,  i.e.  (dP/dd)  =  0.  For  this  limiting  bearing  configuration, 
Reynolds  equation  reduces  to 


dtl  ’  2  3© 1  ’  dz\uM  dzj 


L 


L/D  «  1 
dP/dO  ->  0 


Figure  7:  Short  Length  Bearing  Model. 


(19) 
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The  short  length  bearing  model  provides  (surprisingly)  accurate  results  for  plain  cylindrical  bearings  of 
slenderness  ratios  L/D  <0.50  and  for  small  to  moderate  values  of  the  journal  eccentricity,  e< 0.75  c  [4], 
The  short  length  bearing  model  is  widely  used  for  quick  estimations  of  journal  bearing  static  and  dynamic 
force  performance  characteristics.  Integration  of  equation  (19)  leads  to  the  pressure  distribution 


P(0,z,t)-Pa 


f  ;  Cl) 

e  cos  0  +  e 

<P~- 

sin# 

v  2  J 

C3  H 3 


(20) 


with  H-h/c  -  1  +  s  cos(d)  as  the  dimensionless  film  thickness,  and  e  =  e/c  is  the  journal  eccentricity  ratio; 
[0  <  s  <  1],  £  —  0  means  centered  operation  (typically  a  condition  of  no  load  support),  and  £  —  1.0 
evidences  solid  contact  of  the  journal  with  its  bearing. 

No  lubricant  cavitation  will  occur  if  the  exit  or  discharge  pressure  Pa  is  well  above  the  liquid  cavitation 
pressure.  However,  if  Pa  is  low,  typically  ambient  conditions  at  1  bar,  it  is  almost  certain  that  the  bearing 
will  cavitate  or  show  air  entrainment  when  the  outlet  plenum  is  not  flooded  with  lubricant.  The  cavitation 
model  in  the  short  length  bearing  simply  neglects  any  predicted  negative  pressures  and  equates  them  to 
zero.  This  “chop”  procedure  although  theoretically  not  well  justified  grasps  with  some  accuracy  the  actual 
physics  [5].  Hence,  if  Pa  =  0,  and  from  equation  (21),  the  pressure  field  P>(),  when  cos( 0+a)  <  0.  Thus, 
P>0  in  the  circumferential  region  limited  by 

-<0  +  a<—  -a=0l  <6<e2=  —  -a  (21) 

2  2  2  2 


That  is,  regardless  of  the  type  of  journal  motion,  the  region  of  positive  pressure  has  an  extent  of  n 
(  =  180°);  thus  then  the  infamous  7tfilm  cavitation  model  widely  used  in  the  literature. 


Fluid  film  reaction  forces  on  the  journal  are  evaluated  by  integration  of  the  pressure  field  acting  on  the 
journal  surface.  With  Pa  =  0,  the  radial  and  tangential  forces  (F„  F,)  are  given  by 


/u  RL3 


J 

J\ 


02 


11 

3 

20 

3 


(22) 


where  the  J’s  are  integrals  defined  in  analytical  form  by  Booker  [6].  Note  that  the  fluid  film  forces  are 
proportional  to  the  journal  center  translational  velocities  (e,etp)  as  well  as  the  journal  rotational  speed 

(Q).  The  reaction  forces  depend  linearly  on  the  fluid  viscosity  and  the  bearing  radius  and  grow  rapidly 
with  the  ratio  (L/C)3. 

Hydrodynamic  journal  bearings  are  designed  (and  implemented)  to  support  a  static  load  W,  hereafter 
aligned  with  the  X  axis  for  convenience,  see  Figure  8.  At  the  equilibrium  condition,  denoted  by  a  journal 
center  eccentric  displacement  (e)  with  an  attitude  angle  ( tp).  the  hydrodynamic  bearing  generates  a  reaction 
force  balancing  the  applied  external  load  at  the  rated  rotational  speed  (Cl).  The  equations  of  static 
equilibrium  are 


W  +  Fx  =  0  =>  -W  =  Fx  =  Fr  cos  tp  -Ft  sin  tp 
FY  =  0  =>  0  =  FY  =  Fr  sin  tp  +  Ft  cos  tp 


(23) 
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Static  load 


Figure  8:  Force  Equilibrium  for  Statically  Applied  Load. 


For  static  equilibrium,  e =0,  (f)  =  0,and  0j  —  0  to  02  =  n.  From  equation  (22),  the  static  radial  and  tangential 
film  reaction  forces  are 


//  RL3Q  a2  uRL3  Cl  n-s 

c3  (l-s3)3'  '  4(l -83T 

Figure  9  depicts  the  radial  and  tangential  forces  for  a  typical  short  length  bearing.  The  forces  are 
proportional  to  the  lubricant  viscosity  and  rotor  surface  speed  ( Q.R ),  the  length  (Z/),  and  inversely 
proportional  to  the  radial  clearance  (c2).  Most  importantly,  the  bearing  forces  grow  rapidly  (non-linearly) 
with  the  journal  eccentricity  (s-e/c). 


-  -Fr 

....  Ft 


Figure  9:  Radial  and  Tangential  Forces  for  Short  Length  Bearing. 
p=0.019  Pa.s,  L=0.05  m,  c=0.1  mm,  3,  000  rpm,  L/D=  0.25. 


The  external  load  (W)  is  balanced  by  the  fluid  film  reaction  forces.  Thus, 

Vl6^r  +^-2(l-^2) 


W  =  {F;  +F?)V2  =  /uCIRL 


L 


£ 

4 


(25) 


RTO-EN-AVT-143 


10-11 


Hydrodynamic  Fluid  Film  Bearings  and 

Their  Effect  on  the  Stability  of  Rotating  Machinery 


ORGANIZATION 


and  the  journal  attitude  angle  <f>  is  obtained  from 


tang  (p  = 


Fr 


Xyl(l-£2) 

4  £ 


Note  that  as  the  journal  eccentricity  s  0,  p  —>  n/2,  while  as  s  —>  1,  p  —>  0. 


(26) 


3.1  Design  of  Hydrodynamic  Bearing  -  Selection  of  Operating  Eccentricity 

In  the  design  of  hydrodynamic  journal  bearings,  the  bearing  static  performance  characteristics  are  related 
to  a  unique  dimensionless  parameter  known  as  the  Sommerfeld  Number  (S)  defined  as 


S  = 


//  NLD 


f-Y 


(27) 


where  N  ~(Q/2n)  is  the  rotational  speed  in  revolutions/sec.  In  practice,  the  specific  load  or  pressure  is 
known  as  the  ratio  of  applied  load  to  bearing  projected  area,  i.e.  (W/LD). 


In  short  length  journal  bearings,  a  modified  Sommerfeld  number  (a)  is  defined  and  related  to  ( S )  by  [5,  7]: 


G  =  nS  (L/D)2 


H  D.LR  f  lV 
4W  (,C, 


(28) 


Substitution  of  Eqn.  (28)  into  Eqn.  (25)  relates  the  modified  Sommerfeld  number  to  the  equilibrium 
operating  journal  eccentricity  (e),  i.e. 


iu  Q LRfLY 
4W  U> 


(29) 


At  a  rated  operating  condition,  a  is  known  since  the  bearing  geometry  (R,  L,  c),  rotational  speed  (Q),  fluid 
viscosity  (p)  and  applied  load  (W)  are  known.  Then,  equation  (29)  provides  a  relationship  to  determine 
(iteratively)  the  equilibrium  journal  eccentricity  ratio  s  -(e/c)  required  to  generate  the  fluid  film  reaction 
force  balancing  the  externally  applied  load  W. 

Figures  10  and  1 1  depict  the  modified  Sommerfeld  number  and  attitude  angle  vs.  journal  eccentricity, 
respectively.  Large  Sommerfeld  (a)  numbers;  i.e.  denoting  small  load,  high  speed  Q  or  large  lubricant 
viscosity,  determine  small  operating  journal  eccentricities  or  nearly  centered  operation,  s  —>0,  (j)^>n/2 
(90°).  That  is,  the  journal  eccentricity  vector  is  nearly  orthogonal  to  the  applied  load. 


10-12 


RTO-EN-AVT-143 


Hydrodynamic  Fluid  Film  Bearings  and 
Their  Effect  on  the  Stability  of  Rotating  Machinery 


Sommerfeld  number  for  short  journal  bearing 


journal  eccentricity  (e/c) 


Figure  10:  Modified  Sommerfeld  (?)  Number  versus  Journal  Eccentricity. 


Figure  11:  Equilibrium  Attitude  Angle  versus  Journal  Eccentricity. 


Small  Sommerfeld  (cr)  numbers,  i.e.  denoting  large  load,  low  speed  D.  or  low  lubricant  viscosity, 
determine  large  operating  journal  eccentricities,  s  —>  1.0,  (j)—>0  (0°).  Note  that  the  journal  eccentricity 
vector  is  nearly  parallel  to  the  applied  load. 

Figure  12  shows  the  journal  displacement  within  the  bearing  clearance  for  different  operating  conditions. 
The  journal  eccentricity  approaches  the  clearance  for  large  loads,  low  shaft  speeds  or  light  lubricant 
viscosity,  and  it  is  aligned  with  the  load  vector.  For  small  loads,  high  speeds  or  large  lubricant  viscosities 
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(large  Sommerfeld  numbers),  the  journal  travels  towards  the  bearing  center  and  its  position  is  orthogonal 
to  the  applied  load.  This  peculiar  behavior  is  the  source  of  rotordynamic  instability  as  will  be  shown 
shortly. 


ey/c 

0  0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1 


Figure  12:  Locus  of  Journal  Center  for  Short  Length  Bearing. 


4.0  DYNAMICS  OF  A  RIGID  ROTOR  SUPPORTED  ON  SHORT  LENGTH 
BEARINGS 

Figure  13  depicts  a  symmetric  rigid  rotor  of  mass  2M,  and  supporting  a  static  load  (2Fa-W)  along  the  X 
axis.  The  rotor  is  mounted  on  two  identical  plain  hydrodynamic  journal  bearings.  The  equations  of  motion 
of  the  rotating  system  at  constant  rotational  speed  (Q)  are  given  by  [5]: 


MX  =  Fx  +  M  uCl1  sin(Q  t )  +  Fo 
M  Y  =  FY  +  M  uQr  cos  {Fit) 


(30) 


where  u  is  the  magnitude  of  the  imbalance  vector,  X(t)  and  Y(t)  are  the  coordinates  of  the  rotor  mass 
center,  and  (Fx,  FY) are  the  fluid  film  bearing  reaction  forces.  Since  the  rotor  is  rigid,  the  center  of  mass 
displacements  are  identical  to  those  of  the  journal  bearings,  i.e.  X  (t)  =  ex  ( t ),  Y ( t )  =  eY  ( t ) 
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Static  load 


Figure  13:  Rigid  Rotor  Supported  on  Journal  Bearings, 
(u)  Imbalance,  (e)  Journal  Eccentricity. 


We  are  interested  on  the  rotor  dynamic  behaviour  for  small  amplitude  motions  about  the  equilibrium 
position  defined  by: 


Fx„  -  >  Fr  -  0,  =>  ex  ,  eY  or  e0 ,  (j)c 


(31) 


where  (e0  ,</>„)  denote  the  static  equilibrium  journal  eccentricity  and  attitude  angle, 
respectively.  The  bearing  static  reaction  forces  satisfy 


Fo  +  FXo  =0  =>  F0=FXq=  Frg  cos  (j)0  -  F,o  sin  (f)0 
F yn  =  0  =>  0  =  FYq=  FIq  sin  <f)0  +  Ftn  cos  <f)0 


(32) 


Small  amplitude  journal  motions  about  the  equilibrium  position,  as  represented  in 
Figure  14,  are  defined  as: 


ex  =  ex  +A ex(t),  eY=eY  +  AeY (t) ,  or X  =  X0  +  AX(t),  Y  =  Y0+AY(t) 


(33a) 


or  conversely, 


with 


e(t)  =  e0+  Ae(t),  <f){t)  =  <j)0+  A</>(t) 


dX 
dt 
d2X 
dt 2 


ex  —  A  ex . 

=A«X- 


dY 


dt 

d2Y 


=  eY  =  A  eY 


dt 


2  =«y  =  A«y 


(33b) 


(33c) 
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Figure  14:  Small  Amplitude  Journal  Motions  about  an  Equilibrium  Position. 


The  journal  dynamic  displacements  in  the  (r,  t)  coordinate  system  are  related  to  those  in  the  (X,  Y)  system 
by  the  linear  transformation 


1 

> 

X 

_ i 

COS  (j)0 

-  sin  (j)0 

A  e(t) 

i 

> 

i _ 

sin  (f>0 

COS  <j)0 

1 - 

O 

> 

(34) 


Similar  relationships  hold  for  the  journal  center  velocities  and  accelerations.  Note  that  the  assumption  of 
small  amplitude  motions  requires  A  ex,  A  eY  «c,  i.e.,  the  journal  dynamic  displacements  are  smaller 
than  the  bearing  clearance. 

Recall  that  the  fluid  film  forces  are  general  functions  of  the  journal  center  displacements  and  velocities, 
i.e.  Fa  =  Fa  [ex(t),  eY(t),ex(t),  eY(t)],  a  =  X,Y .  Now,  express  the  bearing  reaction  forces  as  a 
Taylor  Series  expansion  around  the  static  journal  position  (eXo,  eYo),  he. 


Fx  =  FXn 


dX 


dFY 


dY 


dX 


(35) 


Fr  =Fy  +  ^AX  + 

Y  Y°  ex 


dF, 


dY 


Y  A  Y  + 


8F  dF 

— fax+^ay 

dX  dY 


4.1  Definition  of  Dynamic  Force  Coefficients  in  Fluid  Film  Bearings 

Fluid  film  bearing  stiffness  ( K,j)  ij=Xj  and  damping  ( C,j)  y=Xt  Y  force  coefficients  are  defined  as 


jr 

ij  dX , 


r  dF< 

ij  dX 


i,j=X,Y 


(36) 


For  example,  KXY  -  -dFx/dY  coiTesponds  to  a  stiffness  produced  by  a  fluid  force  in  the  X  direction  due  to  a 
journal  static  displacement  in  the  Y  direction.  By  definition,  this  coefficient  is  evaluated  at  the  equilibrium 
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position  with  other  journal  center  displacements  and  velocities  set  to  zero.  The  negative  sign  in  the 
definition  ensures  that  a  positive  magnitude  stiffness  coefficient  corresponds  to  a  restorative  force. 

The  force  coefficients  ( Kxx,  Kyy)  are  known  as  the  direct  stiffness  terms,  while  ( KXy,  Kyx)  are  referred  as 
cross-coupled.  Figure  15  provides  a  pictorial  representation  of  the  bearing  force  coefficients  as  mechanical 
parameters. 


Figure  15:  The  “Physical  Representation”  of  Stiffness  and 
Damping  Coefficients  in  Lubricated  Bearings. 


Inertia  or  added  mass  coefficients  {M,j}ij-Xy  can  also  be  defined  as  M  y=x  y  where  {x ,  Y  r  are 

lj  dXj’ 

journal  center  accelerations.  Inertia  coefficients  are  of  particular  importance  in  super  laminar  and  turbulent 
flow  fluid  film  bearings  and  annular  seals.  The  inertia  force  coefficients  or  apparent  masses  have  a  sound 
physical  interpretation.  These  coefficients  are  always  present  in  a  fluid  film  bearing.  Inertia  coefficients 
can  be  of  large  magnitude,  in  particular  for  dense  liquids.  Flowever,  the  effect  of  inertia  forces  on  the 
dynamic  response  of  rotor-bearing  systems  is  only  of  importance  at  large  excitation  frequencies,  i.e.  high 
squeeze  film  Reynolds  numbers.  (This  fact  also  holds  for  most  mechanical  systems  subjected  to  fast 
transient  motions). 


With  the  given  definitions 

'FAD\ 

Fr(t), 

where  Fx„  =  F0=ViW  and  FYo  =  0.  Note  that  the  defined  force  coefficients  allow  the  representation  of  the 
dynamic  fluid  film  bearing  (or  seal)  forces  in  terms  of  the  fundamental  mechanical  parameters  {K,  C,  and 
M}.  Flowever,  this  does  not  mean  that  these  force  coefficients  must  be  accordance  with  accepted  “physics 
grounded”  knowledge.  For  example,  the  “viscous”  damping  coefficients  may  be  negative,  i.e.  non- 
dissipative,  or  the  stiffness  coefficients  non-restorative. 

From  Eqn.  (30),  the  linear  equations  for  small  amplitude  motions  of  the  rotor-bearing  system  are 


,  the  bearing  reaction  forces  are  represented  as 
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The  literature  presents  the  force  coefficients  in  dimensionless  form  according  to  the  definition: 


7  v  c  n  cQ 

k‘j  =  Kjj  —  ,  Cy  =  Cy  ——  iJ-X,  Y 

^0  ^0 


(39) 


where  Fa  is  the  static  load  applied  on  each  bearing  (along  the  X  direction).  [Note  that  the  total  static  load 
W-2F0is,  shared  by  the  two  bearings  in  a  symmetric  rotor  mount]. 

Lund  [8]  derived  first  the  analytical  formulas  for  the  short  bearing  force  coefficients.  Figures  16  and  17 
depict  the  dimensionless  force  coefficients,  stiffness  and  damping,  as  functions  of  the  journal  eccentricity 
and  of  the  modified  Sommerfeld  number  (a),  respectively.  In  the  figures,  both  representations  are 
necessary  since  at  times  the  journal  eccentricity  is  known  a  priori;  while  most  often,  the  design  parameter, 
i.e.  the  Sommerfeld  number,  is  known  in  advance.  In  general,  the  physical  magnitude  of  the  stiffness  and 
damping  coefficients  increases  rapidly  (nonlinearly)  as  the  journal  eccentricity  increases. 
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Stiffness  coefficients  x  (c/Fo) 
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Figure  16:  Dimensionless  Stiffness  and  Damping  Coefficients 
vs.  Journal  Eccentricity  (e)  for  Short  Journal  Bearing. 
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Stiffness  coefficients  x  (C/Fo) 


-  kyy 
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-  cxx 
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Figure  17:  Dimensionless  Stiffness  and  Damping  Coefficients 
vs.  Sommerfeld  Number  (a)  for  Short  Journal  Bearing. 
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Note  that  the  dimensionless  force  coefficients  do  not  represent  the  actual  physical  trends.  For  example,  at 
e„=0,  Kxx-Kyy-O,  but  the  dimensionless  kxx-kyy  have  non  zero  value.  This  peculiarity  follows  from  the 
definition  of  dimensionless  force  coefficients  using  the  applied  load  (F0).  Recall  that,  as  eo—>0,  the  static 
load  Fa  is  also  naught. 


4.2  Dynamic  Force  Coefficients  for  Journal  Centered  Operation,  i.e.  No  Applied  Load 

As  the  journal  center  approaches  the  bearing  center,  eo—>0,  </>„  -  90°,  and  from  the  formulas  presented, 


K  —K  -k-VnRL3  ” C  -C  ~c~mRL3  K 

/v  XY  YX  *  3  ,  ~  c  ’  ^XX  ^YY  L  3  ~ 

c  4  2  cl 


(40) 


Thus,  at  the  centered  journal  position  a  hydrodynamic  bearing  offers  no  direct  (support)  stiffness  but  only 
cross-coupled  forces.  A  small  load  applied  on  the  bearing  will  cause  a  journal  displacement  in  a  direction 
orthogonal  (transverse)  to  the  load,  as  shown  in  the  schematic  view  below.  This  behaviour  is  common  to 
all  fluid  film  journal  bearings  of  rigid  geometry. 


Non-rotating  structure 


Rotating  structure 


The  significance  of 
the  cross-coupled 
effect  in  fluid  film 
bearings 


5.0  ROTORDYNAMIC  STABILITY  OF  RIGID  ROTOR  SUPPORTED  ON 
SHORT  LENGTH  BEARINGS 

The  linearized  equations  of  motion  are  written  in  dimensionless  form  as  [5,  9] 


"A/ 

cxx 

CXY 

Ax' 

^  XX 

k 

r^XY 

Ax 

2  c* 

cos (r) 

+ 

+ 

=  p  o 

A/ 

_CYX 

CYY  _ 

A  / 

k 

_^YX 

k 

kyy  _ 

_A  y 

sin(r) 

where  Ax=  Ay  =  t  =  Fit,  8  =  —  ,(')  =  — ;  p2  =  M  ^  is  a  dimensionless  mass,  and  kt]  and 

c  c  c  dr’  Fo 

c,,  are  the  dimensionless  stiffness  and  damping  force  coefficients. 

It  is  of  interest  to  determine  if  the  rotor-bearing  system  is  stable  for  small  amplitude  journal  center 
motions  (perturbations)  about  the  equilibrium  position.  To  this  end,  set  the  imbalance  parameter  8  -  0  in 
the  equations  above  to  obtain, 
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k 

^  XX 

l - 

X 

Ax 

"0" 

+ 

k 

_YX 

k 

/Vyy  _ 

Ay. 

0 

(42) 


If  the  rotor-bearing  system  is  to  become  unstable,  this  will  occur  at  a  threshold  speed  of  rotation  (/2t) 
and  the  rotor  will  perform  (undamped)3  orbital  motions  at  a  whirl  frequency  ( cos ).  These  motions, 
satisfying  equation  (42),  are  of  the  form: 

x=Aejt0-'  =Aej”T  ;  y  =B  eiaJ  =  Bej“T;  j  =  V3!  (43) 

where  co  =  coj Qs  is  known  as  the  whirl  frequency  ratio,  i.e.  the  ratio  of  the  rotor  whirl  or  precessional 
frequency  to  the  rotor  onset  speed  of  instability.  Substitution  of  (43)  into  equation  (42)  leads  to  [5]: 

~  Ps®s  +  k-xx  +  j<X>sCx. X  kXy  +  j^sCXY 

kyx  +  j(OsCYX  —  P  s(Os  +  kyy  +  jCO^yy 


pc 

"0" 

|_  B 

0 

(44) 


In  Eqn.  (44),  the  determinant  A,  must  be  zero  for  a  non-trivial  solution  of  the  homogenous  system.  After 
algebraic  manipulation,  the  real  and  imaginary  parts  of  A  render  [5,9] 


2—2  _  t  _  kXxCYY  ~^^YYCX X  ~  C YX^ XY  ~  C XY^YX  _  C  M CO s 
r  s  _  ^  eq  ~  ~  ^ 

CXX  +CJT  Fo 


(45) 


and 


c o :  = 


( keq  ^XX  \^eq  ^ YY  )  ^ XY  '  ^Y 


f  \ 2 

CO „ 


\QSJ 


(46) 


5.1  Threshold  Speed,  Critical  Mass,  Equivalent  Stiffness  and  Whirl  Frequency  Ratio 

For  a  given  value  of  journal  eccentricity  (s0),  i.e.  a  given  Sommerfeld  number  (cr ),  one  evaluates  Eqn. 
(45)  to  obtain  the  equivalent  stiffness  keq,  and  then  Eqn.  (46)  to  get  the  whirl  frequency  ratio  (Os  ■  This 
substitution  then  yields  p2s  =  kn)  j 6) 2  (-critical  mass),  which  in  turn  renders  the  onset  speed  of 
instability  Qs. 

Figures  18  and  19  depict  the  whirl  frequency  ratio  (cc/f2)s  and  the  dimensionless  threshold  speed  of 
instability  (ps)  versus  equilibrium  journal  eccentricity  and  modified  Sommerfeld  number,  respectively. 
Note  that  for  near  centered  journal  operation,  i.e.  large  Sommerfeld  numbers,  the  whirl  frequency  is  0.50, 
i.e.  half-synchronous  whirl. 


Recall  that  in  a  mechanical  system,  an  equivalent  damping  ratio  >  0  causes  the  attenuation  of  motions  induced  by  small 
perturbations  from  an  equilibrium  position.  A  null  damping  ratio  brings  the  system  into  sustained  periodic  motions  without 
decay  or  growth,  thus  denoting  the  threshold  between  stability  and  instability  (amplitude  growing  motions). 
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journal  eccentricity  (ratio) 

high  speed  low  speed 

small  load  <  -  — >  large  load 

large  viscosity  small  viscosity 

Figure  18:  Whirl  Frequency  Ratio  vs.  Sommerfeld  Number  (a)  and  Journal  Eccentricity  (e). 
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Dim  [-]  threshold  speed  (ps) 


Dim  [-]  threshold  speed  (ps) 


high  speed 
small  load  <  - 
large  viscosity 


low  speed 
- >  large  load 

small  viscosity 


Figure  19:  Threshold  Speed  of  Instability  (ps)  vs.  Sommerfeld 
Number  (cr)  and  Journal  Eccentricity  (s). 


On  the  other  hand,  if  one  assumes  that  the  current  rotational  speed  (Q)  is  the  onset  speed  of  instability, 
then  from  the  relations  above  it  follows  the  largest  magnitude  of  Vi  system  mass  (M)  to  make  the  rotor¬ 
bearing  system  unstable.  This  mass  is  known  as  the  critical  mass,  Mc,  and  corresponds  to  the  limit  mass 
which  the  system  can  carry  dynamically.  If  the  rotor  mass  is  equal  to  or  larger  than  twice  Mc,  then  the 
system  will  become  unstable  at  the  rated  speed  Q  4. 

The  whirl  frequency  ratio  (WFR),  coj Qs  ,  is  the  ratio  of  the  rotor  whirl  frequency  to  the  onset  speed  of 
instability.  Note  that  this  ratio,  as  given  in  equation  (46),  depends  only  on  the  fluid  film  bearing 

4  Recall  that  each  bearing  carries  half  the  static  load,  and  also  half  the  dynamic  or  inertia  load  (2.McC  fZ). 
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characteristics  and  the  equilibrium  eccentricity.  The  WFR  is  independent  of  the  rotor  characteristics  (rotor 
mass  and  flexibility)  [5].  Reference  [10]  presents  an  analysis  including  fluid  inertia  effects,  more 
applicable  to  annular  pressure  seals  and  bearings  handling  process  fluids  of  large  density. 

The  parameter  keq  is  a  journal  bearing  (dimensionless)  equivalent  stiffness,  also  depicted  in  Figures  16  and 
17.  From  the  definitions  of  threshold  speed  and  whirl  ratio,  p2s  =  M  Q 2  (C/ Fn )  and  cds  =  <x>s  / Qs  ,  then 


Mco,  =k 


eq 


c 


=  K 


eq 


Thus,  the  whirl  or  precessional  frequency  is 


CO,  = 


=  co„ 


(47) 


(48) 


i.e.,  the  whirl  frequency  equals  the  natural  frequency  of  the  rigid  rotor  supported  on  journal  bearings. 

For  operation  close  to  the  concentric  position,  s„  —>  0,  i.e.  large  Sommerfeld  numbers  (no  load  condition), 
the  force  coefficients  are,  see  equation  (40), 


^■xx  kYy  0,  c  ^  c~yy, 


k  =  -k  ■ 

^XY  YX  ’ 


c  —  c  =  0 

u  XY  ^YX  w 


(49) 


keq  {^XXCXX 


+  cxy^xy)/cxx  —0 


and 


b 

^  XY 
C  XX 


0.50 


as  a  -» 0 


(50) 


This  value  of  whirl  frequency  ratio  (WFR)  is  a  characteristic  of  hydrodynamic  plain  journal  bearings.  The 
WFR  shows  that  at  the  onset  sped  of  instability  the  rotor  whirls  at  its  natural  frequency  equal  to  50%  of  the 
threshold  rotational  speed.  Furthermore,  under  no  externally  applied  loads,  F0= 0,  as  in  vertically 
turbomachinery,  the  bearing  possesses  no  support  stiffness,  i.e.  Keq= 0  and  the  system  natural  frequency 
( co„)  is  zero,  i.e.  the  rotor-bearing  system  must  whirl  at  all  operating  speeds. 

Note  that  if  KXy  -  0,  i.e.  the  bearing  does  not  have  cross-coupled  effects,  then  the  WFR  =  0,  i.e.  no  whirl 
occurs  and  the  system  is  always  dynamically  stable.  Cross-coupled  effects  are  then  responsible  for  the 
instabilities  so  commonly  observed  in  rotors  mounted  on  journal  bearings.  If  the  whirl  frequency  ratio  is 
0.50,  then  the  maximum  rotational  speed  that  the  rotor-bearing  system  can  attain  is  just, 


°mM  =  (T50  =  1(°s  =  2(°n  (51) 

i.e.,  twice  the  natural  frequency  (or  observed  rigid  rotor  critical  speed). 

Figures  18,  19  and  20  depict  the  whirl  frequency  ratio,  the  dimensionless  threshold  speed  (ps)  and  the 
critical  mass  (ps)2  versus  the  Sommerfeld  number  and  equilibrium  journal  eccentricity.  The  results 
demonstrate  that  a  rigid-rotor  supported  on  plain  journal  bearings  is  STABLE  for  journal  eccentricity 
ratios  s  >  0.75  (small  Sommerfeld  numbers )  for  all  L/D  ratios.  Note  that  the  critical  mass  and  whirl 
frequency  ratio  are  nearly  invariant  for  operation  with  journal  eccentricities  (sa)  below  0.50. 
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Critical  mass 


Dimensionless  Critical  mass 


high  speed 
small  load  <— 
large  viscosity 


low  speed 
— >  large  load 

small  viscosity 


Figure  20:  Critical  Mass  (mc=ps2)  vs.  Sommerfeld  Number  (a)  and  Journal  Eccentricity  (s). 


Keep  in  mind  that  increasing  the  rotational  speed  of  the  rotor-bearing  system  determines  larger 
Sommerfeld  numbers,  and  consequently,  operation  at  smaller  journal  eccentricities  for  the  same  applied 
static  load.  Thus,  operation  at  ever  increasing  speeds  will  eventually  lead  to  a  rotor  dynamically  unstable 
system  as  the  analysis  results  show. 

5.2  Effects  of  Rotor  Flexibility  on  Stability  of  System 

A  similar  analysis  can  be  performed  considering  rotor  flexibility  [5,  11].  This  analysis  is  more  laborious 
though  straightforward.  The  analysis  shows  that  rotor  flexibility  does  not  affect  the  whirl  frequency  ratio. 
However,  the  onset  speed  of  instability  is  dramatically  reduced  since  the  natural  frequency  of  the  rotor- 
bearing  system  is  much  lower.  The  relationship  for  the  threshold  speed  of  instability  of  a  flexible  rotor  is: 
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where  the  sub  index  /  denotes  a  flexible  rotor.  Krot  is  the  rotor  stiffness  on  each  side  of  the  mid  disk 
shown  in  the  graph,  and  T  =  Fo  /  K mt  is  the  rotor  static  sag  or  elastic  deformation  at  midspan. 

The  elastic  shaft  and  bearing  are  mounted  in  series,  i.e.  the  bearing  and  shaft  flexibilities  add  (reciprocal  of 
stiffnesses),  and  thus  the  equivalent  system  stiffness  is  lower  than  that  of  the  bearings,  and  therefore  the 
system  natural  frequency  decreases  significantly. 

Figure  21  depicts  the  threshold  speed  of  instability  (psf)  for  a  flexible  rotor  mounted  on  plain  short  length 
journal  bearings.  Note  that  the  more  flexible  the  rotor  is,  the  lower  the  threshold  speed  of  instability.  If  the 
fluid  film  bearings  are  designed  too  stiff  (low  Sommerfeld  numbers),  then  the  natural  frequency  of  the 
rotor-bearing  system  is  just  ( Kro/M)m ,  irrespective  of  the  bearing  configuration. 


Threshold  speed  (ps)  for  flexible  rotor 


T/c=0.1 


•  T/c=l 
-  T/c=10 


low  speed 
large  load  <— 
low  viscosity 


high  speed 
— >  small  load 

large  viscosity 


Figure  21 :  Threshold  Speed  of  Instability  (ps)  for  Flexible  Rotor 
versus  Sommerfeld  Number  (a).  Static  Sag  [Tic)  Varies. 


5.3  Physical  Interpretation  of  Dynamic  Forces  for  Circular  Centered  Whirl 

The  bearing  dynamic  forces  in  the  radial  and  tangential  are 
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Consider  circular  journal  motions  of  amplitude  Ae  at  a  forward  frequency  (oj),  as  shown  in  Figure  22.  At 
the  centered  position,  the  bearing  has  no  direct  stiffnesses,  only  cross-coupled  stiffness  and  direct 
damping,  i.e. 


Krr  =  K„  =  C„  =  C„  =  0 


—  Q  —  —  u RU n 

K  =  Krl=-K,=-C-  C  =  Ctt  =  Crr  =  ~ 


(54) 


Figure  22:  Force  Diagram  for  Circular  Centered  Whirl  Motions. 

And  thus,  the  radial  and  tangential  forces  become 

Frd  =  Ftd  =  -  (Ctt  <*>  -  K rt  )  (55) 

A  destabilizing  force  will  drive  the  journal  in  the  direction  of  the  forward  whirl  motion,  i.e.  Ft> 0  if  the 
equivalent  damping  ( Ceq)  is  negative  (see  Figure  23),  i.e. 

(C„  ~  —  Kn  )  =  C  <0  (56) 

CO 
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whirl 

orbit 


Figure  23:  Forces  Driving  and  Retarding  Rotor  Whirl  Motion. 


At  the  threshold  speed  of  instability,  Kn  =—Cn.  Thus,  unstable  forward  whirl  motions  occur  for  rotor 
speeds  Q  >2cos. 

In  the  (X,  Y)  coordinate  system,  AX-Ae  cos(cot)  and  AY-Ae  sin(cot).  Thus,  the  bearing  dynamic  forces 
become 
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Note  that  (Fx ,  FY)  oppose  the  forward  whirl  motion  for  journal  speeds  fJ  <  2 cos.  For  larger  rotor  speeds, 
the  bearing  forces  become  positive  and  aid  to  the  growth  of  the  forward  whirl  amplitude  of  motion,  as 
shown  graphically  in  Figure  24. 
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„  .  ,  .  .  |  .  . Tr  .  AX<0,  AY>0 

Energy  or  orbit:  Area_orbit  (Kxy-Kyx) 

Figure  24:  Representation  of  Follower  Force  from  Cross-Coupled  Stiffnesses. 


The  work  performed  by  the  bearing  forces  during  a  full  orbital  period  (T-Iti/oj)  equals  [9] 
E  =  | (fXj  AX (t)  +  FYj  A y(f) ) dt  =^FtaAe)dt=-{Cn(o- Krt )coAe2  T  = 

E  =  - -  (2  n  Ae~  )(C„  o)-Krt)=-  2 Areaorbit  Ceq  co 


(58) 


Note  that  E<0  is  equivalent  to  negative  work,  i.e.  energy  removed  or  dissipated  from  the  rotor-bearing 
system.  However  when  E>0,  i.e.  for  Q  >2cos,  the  fluid  film  bearing  adds  "energy"  into  the  rotor-bearing 
system  thus  driving  the  whirl  motion  forward. 

From  this  discussion  one  can  easily  deduce  that  rotor-bearings  evidencing  whirl  orbits  with  skewed  areas 
(sharp  ellipsoids)  will  be  less  prone  to  rotordynamic  instability,  see  Figure  25.  This  type  of  dynamic 
response  is  obtained  by  design  and  construction  of  a  bearing  generating  (direct)  stiffness  asymmetry,  as 
given  in  multiple  pad  bearing  configurations  (elliptic,  multiple-lobe  with  preloads,  pressure-dam 
bearings).  However,  these  bearings  are  limited  to  fixed  orientation  static  loads  and  rotor  spin  in  only  one 
direction. 
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Energy  of  orbit:  Area  x  (Kxy-Kyx) 

Figure  25:  Influence  of  Bearing  Asymmetry  on  Whirl  Orbits. 


5.4  Experimental  Measurements  of  Rotor-Bearing  System  Instability 

The  archival  literature  is  abundant  in  experimental  and  field  descriptions  of  severe  instabilities  induced  by 
fluid  film  bearings  on  rotating  machinery.  As  an  example  of  tests  conducted  at  the  author’s  laboratory  on  a 
high  speed  test  rig,  Figure  26  depicts  recorded  amplitudes  of  motion  versus  shaft  speed  in  a  rigid  rotor 
supported  on  plain  journal  bearings.  The  displacement  measurements  correspond  to  rotor  motions  along 
the  vertical  and  horizontal  planes  (LV,  LF1).  The  curves  with  larger  amplitudes  denote  the  total  amplitudes 
of  motion  while  the  others  in  light  color  show  the  filtered  synchronous  (IX)  motions  with  slow  roll 
compensation.  The  passage  through  a  well-damped  critical  speed  is  evident  at  ~  8.5  krpm.  As  the  shaft 
speed  increases,  the  amplitudes  of  motion  decrease.  Flowever,  at  a  shaft  speed  ~  twice  the  critical  speed, 
the  rotor  becomes  violently  unstable  with  large  amplitude  motions  nearly  equalling  the  bearings’ 
clearances. 


Figure  26:  Amplitudes  of  Rotor  Motion  versus  Shaft  Speed. 
Experimental  Evidence  of  Rotordynamic  Instability. 
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Figure  27  depicts  the  waterfall  of  the  vertical  shaft  motion.  The  graph  shows  the  frequency  content  of  the 
vibration  signal  as  the  rotor  accelerates.  The  synchronous  motions  are  denoted  by  the  IX  line.  The  whirl 
frequency  ratio  is  0.50  at  the  onset  of  the  severe  subsynchronous  motions.  As  the  speed  increases,  the 
whirl  frequency  locks  at  the  system  natural  frequency.  This  phenomenon  is  known  as  oil-whip.  The  rotor 
was  severely  damaged  upon  completion  of  the  experiment. 


Figure  27:  Waterfall  of  Recorded  Rotor  Motion  Demonstrating  Subsynchronous  Whirl. 


6.0  CLOSURE 

Compressors,  turbines,  pumps,  electric  motors,  electric  generators  and  other  rotating  machines  are 
commonly  supported  on  fluid  film  bearings.  In  the  past,  most  applications  implemented  common 
cylindrical  plain  journal  bearings.  As  machines  have  achieved  higher  speeds  and  larger  power,  rotor 
dynamic  instability  problems  such  as  oil  whirl  have  brought  the  need  to  implement  other  bearing 
configurations.  Cutting  axial  grooves  in  the  bearing  to  supply  oil  flow  into  the  lubricated  surfaces 
generates  some  of  these  geometries.  Other  bearing  types  have  various  patterns  of  variable  clearance 
(preload  and  offset)  to  create  a  pad  film  thickness  that  has  strongly  converging  and  diverging  regions,  thus 
generating  a  direct  stiffness  for  operation  even  at  the  journal  centered  position.  Various  other  geometries 
have  evolved  as  well,  such  as  the  tilting  pad  bearing,  which  allows  each  pad  to  pivot,  and  thus  to  take  its 
own  equilibrium  position.  This  feature  usually  results  in  a  strongly  converging  film  region  for  each  loaded 
pad  and  the  near  absence  of  cross-coupled  stiffness  coefficients. 

Tables  2  and  3  summarize  some  of  the  advantages  and  disadvantages  of  various  bearings  in  condensed 
form.  Figure  28  shows  graphical  sketches  for  some  of  the  bearing  configurations  below.  References  [12, 
13,  14]  offer  important  technical  information  on  the  design,  operation  and  stability  considerations  for  the 
most  common  fluid  film  bearings  used  in  industrial  applications,  with  emphasis  in  pumps  and 
compressors. 
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Table  2:  Fixed  Pad  Non-Pre  Loaded  Journal  Bearings 


Bearing  Type 

Advantages 

Disadvantages 

Comments 

Plain  Journal 

1 .  Easy  to  make 

2.  Low  Cost 

1. 

Most  prone  to  oil  whirl 

Round  bearings  are  nearly 
always  “crushed”  to  make 
elliptical  bearings 

Partial  Arc 

1 .  Easy  to  make 

1. 

Poor  vibration  resistance 

Bearing  used  only  on 

2.  Low  Cost 

3.  Low  horsepower  loss 

2. 

Oil  supply  not  easily 
contained 

rather  old  machines 

Axial  Groove 

1 .  Easy  to  make 

2.  Low  Cost 

1. 

Subject  to  oil  whirl 

Round  bearings  are  nearly 
always  “crushed”  to  make 
elliptical  or  multi-lobe 

Floating  Ring 

1 .  Relatively  easy  to  make 

2.  Low  Cost 

1. 

Subject  to  oil  whirl  (two  whirl 
frequencies  from  inner  and 
outer  films  (50%  shaft  speed, 
50%  [shaft  +  ring]  speeds) 

Used  primarily  on  high 
speed  turbochargers  for 

PV  and  CV  engines 

Elliptical 

1 .  Easy  to  make 

2.  Low  Cost 

1. 

Subject  to  oil  whirl  at  high 
speeds 

Probably  most  widely  used 
bearing  at  low  or  moderate 

3.  Good  damping  at  critical 
speeds 

2. 

Load  direction  must  be  known 

rotor  speeds 

Offset  Half 

1 .  Excellent  suppression  of 

1. 

Lair  suppression  of  whirl  at 

High  horizontal  stiffness 

(With 

whirl  at  high  speeds 

moderate  speeds 

and  low  vertical  stiffness  - 

Horizontal  Split) 

2.  Low  Cost 

3.  Easy  to  make 

2. 

Load  direction  must  be  known 

may  become  popular  - 
used  outside  U.S. 

Three  and  Four 

1 .  Good  suppression  of  whirl 

1. 

Expensive  to  make  properly 

Currently  used  by  some 

Lobe 

2.  Overall  good  performance 

3.  Moderate  cost 

2. 

Subject  to  whirl  at  high  speeds 

manufacturers  as  a 
standard  bearing  design 
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Table  3:  Pad  Journal  Bearings  with  Steps,  Dams  or  Pockets,  Tilting  Pad  Bearing 

Bearing  Type 

Advantages 

Disadvantages 

Comments 

Pressure  Dam 
(Single  Dam) 

1.  Good  suppression  of  whirl 

2.  Low  cost 

3.  Good  damping  at  critical 
speeds 

4.  Easy  to  make 

1 .  Goes  unstable  with  little 
warning 

2.  Dam  may  be  subject  to  wear  or 
build  up  over  time 

3.  Load  direction  must  be  known 

Very  popular  in  the 
petrochemical  industry. 

Easy  to  convert  elliptical 
over  to  pressure  dam 

Multi-Dam 

Axial  Groove  or 
Multiple-Lobe 

1 .  Dams  are  relatively  easy  to 
place  in  existing  bearings 

2.  Good  suppression  of  whirl 

3.  Relatively  low  cost 

4.  Good  overall  performance 

1 .  Complex  bearing  requiring 
detailed  analysis 

2.  May  not  suppress  whirl  due  to 
non  bearing  causes 

Used  as  standard  design  by 
some  manufacturers 

Hydrostatic 

1 .  Good  suppression  of  oil 
whirl 

2.  Wide  range  of  design 
parameters 

3.  Moderate  cost 

1 .  Poor  damping  at  critical  speeds 

2.  Requires  careful  design 

3.  Requires  high  pressure 
lubricant  supply 

Generally  high  stiffness 
properties  used  for  high 
precision  rotors 

Non-fixed  Pad  Journal  Bearings 

Bearing  Type 

Advantages 

Disadvantages 

Comments 

Tilting  Pad 
Journal  Bearing 

Flexure  Pivot, 
Tilting  Pad 
Bearing 

1.  Will  not  cause  whirl 
(no  cross  coupling) 

1 .  High  Cost 

2.  Requires  careful  design 

3.  Poor  damping  at  critical 
speeds 

4.  Hard  to  determine  actual 
clearances 

5.  Load  direction  must  be  known 

Widely  used  bearing  to 
stabilize  machines  with 
subsynchronous  non¬ 
bearing  related  excitations 

Foil  Bearing 

1 .  Tolerance  to  misalignment. 

2.  Oil-free 

1 .  High  cost 

2.  Dynamic  performance  not 
well  known  for  heavily  loaded 
machinery 

3.  Prone  to  subsync hronous 
whirl 

Used  mainly  for  low  load 
support  on  high  speed 
machinery  (APU  units) 

Hydrodynamic  Fluid  Film  Bearings  and 
Their  Effect  on  the  Stability  of  Rotating  Machinery 


Tilting  pad  bearing 


LEMON  BORE  BEARING 


Flexure  pivot  hydrostatic  pad  bearing 

Figure  28:  Schematic  Views  of  Various  Radial  Fluid  Film  Bearing  Configurations. 
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